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(§7) Gegenstand der Erfindung ist ein approximates Bildre- 
konstrukttonsverrahren fur Com putertomographiege rata mit 
MehrzeHendete ktoren 1 m Spiral betrieb . Dabei wird zunachst 
fur alie Zeilen <Jes Detektors (2) unabhangig eine Uminterpo- 
lation von Fa char- auf Paralleldaten vorgenomman (sog. 
Rebinning). Mit dtesen Paralleldaten laSt sich eln Bild 
rekonstruieren. indom man jeden Strahl mit einem Gewicht 
zum Bild beitrggen lafct, das vom jeweiligen Abstand des 
Strahls zur Bildebene abhangt Damit erhalt jeder Strahl ein 
Gewicht, das sich im aligemeinen lings des Strahls veran- 
derL indem man diese Gewichtsfunktionen durch eina kurze 
Fourierreihe apprpximiert gelangt man zu einer sehr effizien- 
ten Implementierung des Verfahrens ais Fourierrekonstruk- 
tion* 
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Beschreibung 



Essind Computertomographen bekanrit, bei denen ein Hachendetektor vorgesehen ist, der aus einer Matrix 
von Detektorelementen besteht Er ist demgemaB aus einer Reihe von Zeilen aufgebaut, wobei jede Zefle eine 
5 Reihe von Detektorelementen enthait Der Detektor wird von einem pyramidenf ormigen Rontgenstrahlenbun- 
del getroff en und eriaubt je nach Ausdehnung des Flachendetektors in z-Richtung, d. h. in Richtung der System- 
achse, urn die das MeBsystem aus Rontgenstrahler und Detektor gedreht wird, die simultane Abtastung mehre- 
rer Schichten. Zur Abtastung eines Volumens des Untersuchungsobjektes rotiert das MeBsystem kontinuieriich 
urn die Systemachse, wahrend eine Relativbewegung des Untersuchungsobjektes gegenQber dem MeBsystem in 

io Richtung der Systemachse erfolgt 

Fur die Bildrekonstruktion aus Daten eines Spiralscanners mit einer einzigen Detektorzeile sind in der 
Literatiir verschiedene Verfahren bekannt Beispielsweise wird in der EP 0 430 549 A2 ein Verf ahren angegeben, 
bei dem erne Interpolation in Richtung der z-Achse durchgeffihrt wird, indem man die gewichtete Summe der 
beitragenden Strahlen berechnet, wobei die Gewichte vom Abstand der Strahlen zur Bfldebene abhangen. Bei 

15 einem Mehrzeflendetektor ist der Abstand eines Strahls zur Bildebene wegen der Neigung der Strahlen um 
inren Cone- Winkel nicht mehr fiber den ganzen Strahl konstant, sondern wird zu einer Funktion der betrachte- 
ten Position langs des Strahls. Dadurch wird es erforderlich, die Spiralinterpolation bei der Rfickprojektion fur 
jedes Voxel gestreut durchzufuhren. Die dazu erforderliche Kegelstrahl-Ruckprojektion hat aber einen sehr 
hohen Rechenaufwand. In der DE 44 38 988 Al wird anstelle einer Kegelstrahl-Ruckprojektion eine Gewich- 

20 tung der emzelnen Detektorzeile mit anschlieBender gewohnlicher Rfickprojektion vorgeschiagen. Dabei wird 
jedoch der Cone-Winkel der Strahlen vernachlassigt und es ergibt sich ein Verlust an z-Scharfe und eine 
Zunahme von Cone-Beam- Artefakten. 

Das erfindungsgemaBe Verfahren eriaubt die Berficksichtigung der Neigung der Strahlen (Cone-Winkell 
ohnedaB eine Kegelstrahl-Ruckprojektion erforderlich ist 

25 Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zur Bildrekonstruktion bei einem Mehrzeilendetek- 
torrComputertomographen im Spiraibetrieb anzugeben, das eine genaue BUdrekonstruktion innerhalb kurzer 
Zeit ermSglicht. 

Diese Aufgabe ist erfindungsgemaB gelost durch die Merkmale des Patentanspruches. Das erfindungsgemaBe 
Verfahren ist em approximatives Verfahren und stellt eine sehr effiziente Implementierung als Fourier-Rekon- 
30 struktion dar. Dieses Vorgehen vermeidet die nicht korrekte Vertauschung von z-Interpolation und Faltung in 
Gleichung 21 und eriaubt die Berficksichtigung der Neigungswinkel der Strahlen, ohne daB eine aufwendige 
Kegelstrahl-Rfickprojektion erforderlich ware. 

Die Erfindung ist nachf olgend anhand der Zeichnung naher erlautert Es zeigen 

Rg. 1 die wesentlichen Teiie eines Computertbmographen zur Eriauterung des Erfindungsgedankens und 
Fig. 2 bis 4 geometrische Darsteflungen zur Eriauterung des Bildrekonstruktionsverfahrens bei dem Compu- 
tertomographen gemaB Fig. 1. 

In Fig. 1 ist mit 1 der Fokus eines Rontgenstrahlers 9 und mit 2 der Detektor eines Comptrtertomographen 
bezeichnet. Der Detektor 2 besteht aus einer Reihe paralleler Zeflen, von denen jede eine Reihe von Detektor- 
elementen enthait Eine Primarstrahlenblende blendet ein pyramidenformiges Rdntgenstrahlenbfindel ein, das 
genau auf dem Detektor 2 auftrifft Von diesem Rontgenstrahlenbfindel sind in der Fig. 1 beispielhaft nur der 
Mittelstrahl 4 und ein weiterer Strahl 3 dargesteUt Zur Abtastung eines in einem MeBfeld 5 Iiegenden Untersu- 
chungsobjektes rotiert das MeBsystem aus Rontgenstrahler 9 und Detektor 2 um eine Systemachse a Die dabei 
gebildeten Detektorsignale werden einem Rechner 7 zugefuhrt, der ein Bfld des untersuchten Volumens des 
Patienten berechnet, das auf einem Monitor 8 wiedergegeben wird. Fur die Volumenabtastung erfolgt eine 
45 lattvbewegung in Richtung der Systemachse 6 (z-Richtung) zwischen dem MeBsystem 1, 2, 9 und dem im 
MeBfeld 5 Iiegenden Untersuchungsobjekt 



Definition der Geometrie 



50 u g v! erfautert schematisch die Geometrie des zugrundegelegten Mehrzeilendetektor-Computertomogra- 
phen. Der Detektor 2 ist als Zylinderoberflache ausgebfldet Der Radius dieses Zyhnders ist Ri + R<j. Der Fokus 
t liegt auf der Zylinderachse. Der Detektor 2 und der Strahler 9 sind starr miteinander verbunden. Wahrend der 
Aufhahme bewegt sich der Strahler 9 relativ zum zu untersuchenden Objekt auf einer Spiralbahn mit Radius Rf. 
Der Abstand des Detektors 2 vom Drehzentrum sei R4 

55 Der Projektionswinkel ist mit a bezeichnet Innerhalb einer Projektion ist ein Strahl durch seine Parameter (J 
und 5oet eindeutig bestimmt Dabei ist p der Facherwinkei des Strahls und Q>ct die z-Position des zu ihm 
gehorenden Detektorelements relativ zur z-Position des Fokus. 



&>et *= ZDct - ZF(a) (1) 

Die spiralformige Fokusbahn hat die Steigung s und ist gegeben durch 
zp(a) = zf* + a • s (2) 

65 zf,o ist die z-Position der Spirale fur a = 0. 

Da naturiich vom System nur diskrete MeBwerte aufgenommen werden konnen, sind a, p und tp* diskret: 



2 
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a n = nAa, mit Aa = tt— und b = 0...(WL-1), (3) 

^2* * 

P m = (m + 0,5+AJkOAp, mit m = ~f —(f ~ ^ * und (4) 
q--Zp + 0J5+AQ)AZ, Det , mit ? = Om^-I) . (5) 



Der bei an, p m und £o e t^ gemessene Strahl tragt den MeBwert p(ctn, Pm, Coet. q) oder kurz p(n, m, q). 

Fig. 2 zeigt die Projektion der Anordnung aus Fig. 1 in die x-y-Ebene. Die beiden Variablen a und P kenn- 
zeicfanen eindeutig die Projektion eines Strahls in diese Ebene. Den gleichen Zweck erfQUen auch die beiden 
Variablen 0 und p, die mit a und p fiber die beiden folgenden Gleichungen verknupft sind: 

9 = a + P + | # und (6) 

p = -Rfsinp (7) 

Wahrend a und P zur Beschreibung der Strahlen in Fachergeometrie verwendet werden, dienen die Parameter 
9 und p der Beschreibung der Strahlen in ParaUelgeometrie, nach der Durchfuhrung des im folgenden erklarten 
Rebinning, in dem von Facher- auf eine spezielle ParaUelgeometrie uminterpoliert wird. 

Das Rebinning - 



tf ar ft?** \ 
p k = (k + 0£+AK)Ap , mit k = ~— ~( — -M* ^id 



(8) 



e^/Ae, mit i = o...(iv5T r -i) oj 



10 



Aa ist das Projektionswinkelinkrement zwischen aufeinanderfolgenden Projektionen. Np^ ist die Zahl der 
wahrend eines Vollumlaufs aufgenommenen Projektionen, Np ist die Zahl der insgesamt vorliegenden Projek- 
donen. N ist die Zahl der Kanale pro Zeile. Der Einf achheit halber beschranken wir die folgenden Betrachtungen 15 
auf den Fall geradzahliger N. AM ist das sog. Alignment im FacherwinkeL 

A^Det ist der Abstand zweier Detektorzeilen in z-Richtung, Nrows ist die Zahl der Detektorzeilen und AQ ist 
ein mogliches Alignment in der z-Richtung. Die in das Drehzentrum skalierte ^Coordinate £o e t wird mit £ 
bezeichnet Mithin gilt. 
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In diesem ersten Verarbeitungsschritt des Algorithmic werden die Mehrzeilen-Facherdaten in ParaUeldaten 
umsortiert und uminterpoliert Dieser Schritt wird fur alle Detektorzeilen unabhangig voneinander durchge- 
fuhrt Dabei wird zunachst die Schraglage der Strahlen sowie ihre jeweils verschiedene z-Position vernachlas- 45 
sigt, & h. die Strahlen werden zunachst so behandelt als seien sie in einem zweidimensionalen CT-Gerat in einer 
Ebene z = const gemessen worden. Erst spater wird anhand der genauen geometrischen Lage und Orientierung 
der Originalstrahlen die den erzeugten ParaUelstrahlen zuzuordnende Lage und Orientierung ermittelt 

Das Rebinning ist nachfolgend anhand der Fig. 2 genauer erklart Ein Strahl in ParaUelgeometrie ist gegeben 
durch 6 und p. In den ParaUelprojekdonen Uegen Abtastwerte an den diskreten SteUen 50 



55 



vor, wobei ubUcherweise N*" « N und N p par = Np gesetzt wird, d. h. man berechnet die gleiche Zahl an 
Parallelprojektionen mh der gleichen Anzahl Kanale wie in Fachergeometrie. 60 

Der Begriff Paraflelprojektion bezeichnet hier einen Satz paraUeler Strahlen mit gleichem 9i und q. 

Das Rebinning ist im wesentlichen beschrieben durch (6) und (7). Auflosen dieser Gleichungen tiefert eine 
Vorschrift fur die Auswahl der Facherstrahlen, die gesuchten Parallelstrahlen entsprechen. Man erhalt: 
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und (10) 
& = 0,-2-0. (ID 

AuBerdem wird gesetzt: 
q = const (12) 

Die gewunschten Abtastwerte in den Parallelprojektionen konnen anhand der Gleichnngen (10), (11) und (12) 
errechnet werden. '- 
15 Die Berechnung eines Paralleldatums p(©i, pk, p) oder kurz p(I, k, q) aus den Facherdaten p(n, m, q) wird wie 
folgtvollzogen: 

Zunachs t werden aus Gleichung 1 0 und Gleichung 1 1 die Iadizes ft, in und q des Facherstrahls ermittelt, der dem 
Strahi 1, k, q in ParaHelgeometrie entspricbt Es ergibt sich: 



Aa 
und. 

A 

AP 



sowie 
A = q 



35 Daft und m im allgemeinen nicht ganzzahlig sind, muB zwischen den nachsten Nachbarn interpoliert werden. 
Fur den Spezialf all bilineare Interpolation erhalt man mit 

nio = floor (ft), nw *- ceil (ft) 
mio = floor (ft), mm = ceil (ft) 



pdtq) =» (nu — ft)[(mhi - m)p(nio, mu>, q) + (m - mio)p(nio, mhi,q) + (n - nio)[(mu - ro)p(nhi, mio, q) + (m - 
mi 0 )p(nhi v mhb q)] 



Im f olgenden wird die Lage und Orientierung der erzeugten Paraflelstrahlen besdmmt 
45 Das p-t-Koordinatensystem in Fig* 2 sei gegenuber dem x-y-System urn den Wlnkel urn die z-Achse 
verdrehL p bezeichnet die Position eines Strahls innerhalb einer Parallelprojektion. t bezeichnet eine Position 
auf einem Strahi, d h. t ist eine {Coordinate in Strahllangsrichtung. Die Stelle t = 0 entspricht dem Schnittpunkt 
des Strahls mit einer Geraden, die sowohl auf dem Strahi als auch auf der z-Achse senkrecht steht, <L h. der Punkt 
t = 0 ist derjenige Punkt auf dem Strahi mit dem kieinsten Abstand zur z-Achse. Die z-Position eines Strahls bei 
50 t = 0 ist ein weiterer die Lage des Strahls beschreibender Parameter. Wir bezeichnen ihn mit z t *=o. Er ist 
gegebendurch 



55 z t=o(fii,Pk,q)_ - + C |l~ . (13) 



4 



Dabei ist zf(&i, pk) die Fokus z-Position der Strahlen (9u pkX Man kann sie errechnen, indem man zunachst (10) 
60 in(ll)einsetzt: 




(14) 

und dann (2) anwendet Das Ergebnis ist: 
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*F&l>Pd - 0 • (0 f - 1 + asin^)) /. (15) 



2 

' ) 

Wegen des asin(.)-Terms 1st die z-Position der Strahlen innerhalb einer Parallelprojektion kanalabharigig. 
Allerdings haben alle Strahlen einer Parallelprojektion den gleichen Cone-winke^ (das ist der Winkel der 
Strahlen gegenuber der x-y-Ebene). 

Aus der Vielzahl der erzeugten Parallelprojektionen wird nun im folgenden ein Bild rekonstruiert 

Definition der Rekonstruktionsvorschrift 



(v) = i 1 ^Serade v 
[ -1 fur ungerade v 

Um die Notation zu vereinfachen, fuhren wir modifizierte Projektionen und Gewichtungsfunktionen ein, 
denen fur ungeradzahlige v die Koordinatenrichtungen umgekehrt sind. 

P* Zimg* l & V- q) = K(r(y)t, r(v)p, r^, Z 0 , V, q) , ( 18 ) 

p(6 lo + vie, p, q} = + Vic, r(v)p, • ( 19 > 

Dabei ist der Abstand des Strahis in z-Richtung von der Bildebene bei zbng gegeben durch 



m 



s 



10 



Pie Projektionen aller Strahlen der Parallelprojektionen zu den Winkeln 9 = 9*,+ vn (v beliebig ganzzahlig) 
in die Bildebene sind parallel, wobei die Strahlrichtung fur geradzahlige v entgegengesetzt zu der fur ungerad- 
zahlige v ist Bei der Berechnung des BOdes summieren wir daher nur fiber einen Halbumlauf (Projektionswin- is 
kelbereich 



20 



berucksichtigen aber fur alle lo jeweils alle Parallelprojektionen mit 9 = % + vrc fur alle v und q. Abb. 3 
veranschaulicht ein Abtastmuster eines Mehrzeilendetektors als Schnitt entlang der z-Achse, d h. die Zeichen- 
ebene ist die z-t-Ebene. Die dargestellten Strahlen sind die bei pk =» 0. 

Ahnlich wie bei vielen anderen approximativen Verfahren soil in unserer Definition des zu rekonstruierenden 25 
Bildes das Gewicht, mit dera ein Strahl zum Bfld beitragt, vom Abstand in z-Richtung des Strahis zur Bildebene 
abhangen. Dieser Abstand verandert sich jedoch im allgemeinen entlang eines Strahis, so daB das Gewicht, mit 
dem ein Strahl zur Bildebene beitragt ebenso in Strahllangsrichtung variiert Indem man fur ein festes lo, alle 
beitragenden Parallelprojektionen p(9fe + vjt, pi& q) (fur alle (v, q)) ink ihrea zugehorigen Gewichtsfunktionen 
versieht und aufaddiert erhalt man ein gewichtetes Mittel aus diesen Parallelprojektionen mit sich sowohl in 30 
Strahllangsrichtung als auch in Kanalrichtung verandernden Gewichten. Diese gemittelte Parallelprojektion 
wird als Approximation einer unter dem Projektionswinkel 6io vollstandig in der. Bildebene gemessenen Parallel- 
projektion verwendet Die Gewichtungsfunktionen mfissen aus Nonnierungsgrunden so beschaffen sein, daB 
sich die Gewichte aller aufzusummierenden Strahlen jeweils zu 1 erganzen. Die auf eine bestimmte Parallelpro- 
jektion anzuwendende normierte Gewichtsfunktion h n (.) ergibt sich, wenn man die abstandsabhangige Ge- 35 
wichtsfunktion hQ durch die Summe aller beitragenden Gewichte teilt: 

K&p,z imr l Q *V,q) = = * — * . (16) 40 

£ft(^(r(v)-ftr(v)^,^,^v, ? )) 

Weil die Orientierung der Parallelprojektionen fur gerade bzw. ungerade v gerade entgegengesetzt ist, 
muss en fur eine der beiden Falle die beiden Koordinatenrichtungen p und t umgedreht werden. Dies wird durch 45 
r(v) bewerkstelligt Es gilt: 
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Unter Benutzung dieser gemittelten Parallelprojektion kann nun ein gewohnliches Parallelrfickprojektions 
bilddefiniert werden: ■ „ 



5 /fcy>--7^X2Jl^ft*^^ v 'ff^ +v *^« ) | - (21) 



10 



L(- jccos 0^ - y sin 6^ -p fc ) 

15 

Der Term in eckigen Klammern ist die gemittelte Parallelprojektion. 

Die Abhangigkeit von t in ESQ berucksichtigt die Schraglage der Strahlen. Die Auswertung der vorigen 
Gleichung wurde folgende Schritte erf ordem: 

20 — Gewichtete Summation fiber v und q fur alle t 

— Faltung der gewichteten Summe mit L(.) 

— Zuschlagzum bet rachteten Voxel x,y,z- 

Da sich die t- und die p- Abhangigkeit in Ed(.) nicht separieren lassen, kann die Reihenfolge von Summation 
25 uber (v, q) und Faltung fiber k nicht vertauscht werden. Es muS also strenggenommen fur jedes t und damit fur 
jedes Pixel eine eigene Faltung durchgefuhrt werden, Der Rechenaufwand dafur ist aber fur eine praktische 
Anwendung zu hoch. Ffihrt man die Faltung dennoch vor der gewichteten Summation aus, so kann die Interpola- 
tion in z-Richtung bei der Ruckprojektion erfolgen. Man begeht dazu allerdings einen Fehler, Das erfindungsge- 
maBe Verfahren vermeidet das unzulassige Vertauschen der Reihenfolge von gewichteter Summation und 
30 Faltung ohne einen hohen Rechenaufwand zur Folge zu haben. Das Verfahren kommt ohne die aufwendige 
Kegelstrahl-Ruckprojektion aus und berficksichtigt dennoch die Cone-Winkel der Strahlen. 

Durch periodische Wiederholung der gefalteten gemittelten Projektion erreicht man, daB das Spektrum in 
Kanalrichtung diskret wird. Das IJild ist im interessierenden Bereich von dieser periodischen Wiederholung nicht 
betroffen, wenn man davon ausgeht, daB die gefalteten Projektionen endliche Ausdehnung haben und die 
35 Periodewentsprechendgewahlt wird Das Ergebnis ist: 



40 



45 



65 



Innerhalb einer Parallelprojektion ist h^.) eine fur jeden Kanal pk verschiedene eindimensionale Funktion von 
50 t Fur jede Bildebene zim g und fur jeden Wert v, q und pk erhalten wir also eine unterschiedliche Funktion Efo 
. . .), die wir im voraus berechnen und abspeichern konnen. Start eine TabeUe der Funktioiiswerte abzuspeichem, 
approximieren wir die Gewichts funktion mit einer Fourierreihe und speichern die Fourierkoeffizienten. Wegen 
der speziellen Symmetrieeigenschaften der spiralformigen Fokusbahn muB diese TabeUe nicht fur verschiedene 
Projektionswinkei abgespeichert werden. Eine Anderung des Projektionswinkels entspricht namlich einer An- 
55 derung von Zimg. Im folgenden ist erlautert, wie (22) im Frequenzbereich implementiert werden kann, was zu 
einer sog. Fourierrekonstruktion ffihrt 

Implementierung im Frequenzbereich 

60 Zunachst approximieren wir die Gewichtungsfunktionen wie oben erwahnt durch Fourierreihen: 

— J, i^a 



Dadurch werden die Gewichtsfunktionen p eriodisch. Die Periode £J muB geeignet gewahlt werden. 
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Setzt man dies in (22) ein, so erhalt man zunachst: 

k n * 

2*Z< c vSP* Zimr Z 0» v > 9 )e p{\ + Vic, Pk , q) 
£(- xcosG^ - y aiie^ -p k - i\w) 

In dieser Form laBt sich leicht das Spektrum des Bfldes durch zweidimensionale Fouriertransformation is. 
berechnen. Man erhalt: 

2/w 

(25) 20 



_ 2j 
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Die Auswertung dieser Gleichung beinhaltet folgende Schritte: Multiplikation der Projektionen p(9i„ + vj^ pk, 35 
q) mitd en Gewichtungskoeffizienten Cn(pk, Zung, lo v, q) und Summation uber alle beitragenden (v, q). 

FFT der aufsummierten, mit den Koeffizienten muitiplizierten Projektionen entlang der pk-Richtung, und 
zwar fur aHe u, 

Multiplikation des so erzeugten Spektrums mit LQ, wiederum fur alle u. 

Die obengenannten Schritte sind fur alle lo durchzufuhren. 40 
Abb. 4 zeigt das Spektrum einer einzelnen gemittelten Projektion.unter einem bestimmten ProjektiorewinkeL 

Das Spektrum des Bildes ergibt sich durch Oberlagerung aller dieser Behrage uber einen Halbumlaut 
Das resultierende Spektrum kann durch verschiedene Techniken in ein kartesisches Raster uminterpoliert 

(z. B. H. Schomberg, J. Hmmer: The Gridding Method for Image Reconstruction by Fourier Transformation 9 , 

IEEE Tr ansa ctions on Medical Imaging, V6L14, No3, September 1995) und dann mitt els zweidimensionaler 45 

inverser FFT in den Ortsraum zurQcktransformiert werden. 

Patentanspruch 

Approximatives Verfahren zur BOdrekonstruktion bei einem Mehrzeilendetektor-Computertomographen 50 
im SpLralbetrieb, bei dem die vom Detektor (2) gelieferten Daten zunachst fur jede Detektorzeile getrennt 
von Facher- auf Parallelgeometrie umsortiert und uminterpoliert werden, anschlieOend ein Bild aus der 
Menge der erzeugten Parallelprojektionen berechnet wird, wobei das Bild so definiert ist, daB jeder Strahl 
(3, 4) mit einem Gewicht zum Bild beitragt, das von seinem jeweOigen, sich in StrahUangsrichtung indernden 
Abstand zur Bildebene abhangt, der durch diese Abstandsabhangigkeit entstehende Gewichtsverlauf jedes 55 
Strahls (3, 4) in StrahUangsrichtung schlieQlich durch erne FourierreOxe mit den Koeffizienten ji = — Nji 
, Np, approximiert wird, tndem zunachst jeder Kanal einer jeden beitragenden Parallelprojektion mit den 
zugehdrigen qi multipliziert wird, sodann uber alle beitragenden Parallelprojektionen summiert und an- 
schlieBend fur alle p. eine FFT in Kanalrichtung durchgefuhrt und das Ergebnis fur aHe \i mit dem Spektrum 
des Faitungskerns multipliziert wird, wobei das sich ergebende Spektrum der zu einem Projektionswinkel 60 
gehdrenden Bildkomponente damit nur Werte auf diskreten Stutzstellen hat, die auf 2Nn + l paraflelen 
Linien liegen, und das Gesamtspektrum sich durch Summation uber alle Projektionswinkel ergibt und 
schlieBlich durch bekannte Gridding-Techniken auf ein kartesisches Frequenzraster uminterpoliert und 
anschliefiend in den Ortsbereich (Bild) zurucktransf ormiert wird 
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